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Preface 


When  trying  to  make  a numerical  forecasting  model  for  the  oceanic 
circulation,  the  lack  of  knowledgeable  boundary  conditions  at  the  lateral 
boundaries  of  an  arbitrary  body  of  water  was  acutely  felt.  Parts  I-  V of 
this  report  is  an  attempt  to  remedy  that  for  the  North  Atlantic.  This, 
however,  suffers  from  many  defects:  The  model  is  not  all  that  one  de- 
sires; the  resolution  is  low  and  the  equatorial  region  is  not  properly 
treated.  It  is,  however  hesitatingly,  introduced  as  a first  attempt  to  be 
used,  and  abused  and  eventually  to  help  construct  a better  model. 

Part  VI  deals  with  two  side  problems  that  were  of  interest  during 
this  work,  and  were  judged  of  general  interest  to  warrant  reporting  them. 

The  most  interesting  part,  will  probably  be  Part  VII  which  contains 
the  prediction  model  and  tests  based  on  actual  observations.  This  will 
appear  later  in  1964. 

The  U.  S.  Naval  Oceanographic  Office  has  generously  and  sympa- 
thetically sponsored  this  work  under  Contract  N62306-794. 

The  staff  of  the  AEC  Computing  Center  at  New  York  University 
were  very  helpful  and  understanding. 

To  Mrs.  Lillian  Bloom  goes  the  thanks  of  the  authors  for  her  neat 
typing  of  a difficult  piece  of  work  written  almost  in  hieroglyphic.  The 
authors  also  wish  to  thank  Miss  Roberta  Bloom,  Mrs.  Gertrude  Fisher 
and  Mr.  Soliman  Lotaief  for  their  help  in  preparing  the  figures  which  ap- 
pear in  Part  VI  of  this  report. 


Abstract : 


Starting  with  the  wind  stress  over  the  North  Atlantic  Ocean  com- 
puted by  Hidaka  (1958),  a model  of  the  three  dimensional  picture  of  the 
horizontal  velocity  is  developed.  This  is  presented  as  a possible  cli- 
matological atlas  for  the  currents.  Limitations  of  the  method  are  dis- 
cussed. 

Introduction: 

While  attempting  to  forecast  the  oceanic  circulation,  the  problem 
of  lateral  boundary  conditions  had  to  be  tackled.  Any  oceanic  area  not 
involving  the  whole  world  ocean,  will  have  at  least  one  lateral  water 
boundary.  Specifying  the  exchange  of  water  across  this  boundary  affects 
the  circulation  inside  the  area  basically.  In  general,  to  know  the  dis- 
tribution of  the  water  exchange  across  any  arbitrary  water  boundary  in 
a body  of  water,  is  equivalent  to  knowing  the  circulation  everywhere  in 
that  body.  If  this  is  known  for  all  time,  this  implies  that  the  prediction 
of  the  circulation  is  known.  The  circle  of  argument  is  thus  complete: 

To  make  a prediction  of  the  circulation  in  an  arbitrary  are.  in  the  ocean 
the  boundary  exchange  should  be  known,  and  to  know  the  boundary  ex- 
change at  all  time  means  that  a prediction  must  have  been  made.  As 
a compromise  solution,  oceanic  circulation  was  treated  as  a steady  state 
problem,  the  water  being  driven  by  wind  stress.  The  North  Atlantic  was 
chosen  to  compute  numerically  the  currents  using  the  wind  stress  pub- 
lished by  Hidaka  (1958).  These  were  published  by  season,  and  for  the 


annual  mean,  and  correspondingly,  five  circulation  patterns  were  cal- 
culated . 


Fundamental  Equations  and  Assumptions : 

The  equations  used  were  of  the  form: 
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and  u - v - 0 at  the  lateral  boundaries  and  at  the  bottom. 
The  meanings  of  the  symbols  are  as  follows: 


0 Rate  of  change  of  Coriolis  parameter  with  latitude. 
^ North  latitude. 

X East  longitude, 

p Water  density. 

T Wind  stress  in  c.g.s.  units. 

Tx  Zonal  component  of  wind  stress. 

T Meridional  component  of  wind  stress. 


O)  Angular  velocity  of  the  earth's  rotation, 
f Coriolis  parameter  = 2 to  sin  <p  . 
h Constant  depth  of  the  ocean, 
p Pressure, 
u Zonal  velocity, 
v Meridional  velocity, 
x Zonal  coordinate  - R X cos  <p . 
y Meridional  coordinate  = R <p. 
z Vertical  coordinate  increasing  downwards. 

Au  Exchange  coefficient  in  the  horizontal  direction. 

Ay  Exchange  coefficient  in  the  vertical  direction. 

R The  radius  of  the  earth. 

Conditions  imposed  by  the  equations  are: 

(1)  The  flow  is  steady. 

(2)  Vertical  velocity  is  equal  to  zero. 

(3)  Nonlinear  terms  are  negligible,  if  compared  with  terms  re- 
tained in  the  equations. 

(4)  The  exchange  coefficients  are  constant. 

(5)  The  fluid  is  incompressible  and  the  flow  is  nondivergent. 

Setting  the  equations  for  solution: 

When  (1)  is  differentiated  partially  with  respect  to  y and  (2)  with 
respect  to  x,  and  the  first  differentiated  equation  subtracted  from  the 
second,  then  equation  (5)  result: 
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The  velocity  components  u and  v are  then  expanded  in  Fourier  series 
such  that: 


u - £ u.  cos  ~ ^ 
i = 1 J 2h 


v = 2 v.  cos  ^ ' 1^ffz 
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This  expansion  automatical  satisfies  the  boundary  condition  at  the  bottom, 
and  can  represent  any  natural  profile  in  the  ocean. 

Then  it  can  be  shown  from  the  boundary  conditions  at  the  surface 
and  bottom  that: 
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When  substitutions  from  (1),  (6)  and  (7)  are  made  i 
results 


into  (5),  equation  (8) 
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As  all  dependency  on  z in  (8)  appears  in  the  term  cos  and  as 

2h 

h ) is  true  for  all  z,  it  follows  that  the  coefficients  of  the  cosine  term 
can  be  equated  to  zero  independently.  Thus,  equation  (8)  can  be  replaced 
by  the  system  (8'). 


From  (3)  a stream  function 

J 


can  be  introduced  such  that: 
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Each  equation  of  the  system  (8')  can  now  be  set  separately  for  solution. 
Thus,  dropping  j when  it  appears  as  a subscript,  (9)  represents  a typical 
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equation  of  the  system  (8'). 
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Equation  (9)  can  be  put  in  the  finite  difference  form  by  using  the  follow- 
ing approximations  consistent  with  the  retention  of  the  4th  order  deriva- 
tives 
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Computations  and  Results: 

The  values  used  for  the  parameters  were  as  follows: 


Av 

= 

c n -1  -j 

50  gm  cm  sec 

ah 

= 

in8  -1 

10  gm  cm  sec 

h 

= 

2000  meters 

Ay 

= 

555.6  km 

Ax 

= 

Ay  cos  <p 

R 

= 

6371  km 

P 

= 

1.027  gm  cm  8 

60 

— 

7.  292  X 10"5  sec' 

Values  for  wind  stress  were  taken  from  the  computations  by 
Hidaka  (1958)  these  values  appear  at  the  beginning  of  the  different  parts 
of  this  report.  Part  I contains  the  annual  mean,  and  the  seasons  start- 
ing by  the  spring  appear  in  Parts  II  through  V.  The  wind  stress  curl, 
defined  as  the  right  hand  side  of  equation  (10)  is  then  given.  Solving  the 
system  of  equations  (10)  as  a system  of  simultaneous  linear  algebraic 
equations  was  found  to  be  more  accurate  and  less  time  consuming  than 
solving  it  by  iteration.  The  number  of  grid  points  was  147  or  less,  and 
the  resulting  matrix  was  inverted  in  less  than  one  minute.  The  curl  of 
the  wind  stress  could  not  be  accurately  determined  near  the  boundaries 
because  values  on  land  were  taken  as  zero.  The  first  hundred  Fourier 
coefficients  are  given  for  the  grid  points,  and  stream  functions,  zonal 
and  meridional  velocities  for  selected  levels  based  on  the  hundred  Fourier 
coefficients  are  then  given.  The  depth  of  the  selected  levels  appear  in 
Table  I.  The  highest  velocities  appeared  at  the  surface,  but  secondary 
maxima  occurred  at  subsurface  depths.  The  greatest  velocity  in  any 
season  did  not  exceed  100  cm/sec,  and  occurred  near  the  western  bound- 
ary. This  is  less  than  half  the  recognized  value  for  the  western  current, 
but  remembering  that  the  grid  size  is  5 ‘square,  so  that  this  value  is  an 
averaged  value,  this  value  seems  reasonable. 

Discussion  ; 

The  model  discussed  here  is  proposed  as  a probable  climatic 
picture  of  the  North  Atlantic  circulation.  It  should  be  tested  by  compar- 
ing it  to  observations,  and  if  found  adequate,  it  can  serve  as  a skeleton 


9 


to  help  interpret  the  rapidly  accumulating  amount  of  observations.  It  ex- 
hibits the  main  observed  surface  features,  e.g.  the  northern  current 
system,  the  Gulf  Stream  system,  and  the  equatorial  current  system.  Cer- 
tain differences  are  observed,  however,  like  the  north  going  current  off 
the  African  coast.  The  subsurface  currents  cannot  be  verified  because 
of  lack  of  observations,  but  the  under  current  below  the  Gulf  Stream, 
which  was  observed  sometimes,  is  clearly  shown. 

Modifications  of  the  profiles  in  the  model  can  be  accomplished  by 
changing  the  parameters  Ay  and  A^.  Computations  have  been  repeated 
taking  for  Ay  and  A^  the  pairs  of  values  100,  10  ; 100,  10  ; 500,  10  . 

No  significant  difference  in  the  pattern  occurred,  but  differences  occurred 
in  details  e.g.  The  position  and  magnitude  of  the  maximum  stream  func- 
tion at  different  depths.  The  use  of  exchange  coefficients  varying  in  space 
was  not  done,  as  precise  knowledge  of  their  magnitude  is  lacking. 

The  greatest  restriction  of  the  model,  however,  is  the  forced  ab- 
sence of  vertical  velocity.  In  this  model,  it  is  a price  paid  to  enable  a 
three  dimensional  picture  of  the  horizontal  velocities.  It  is  only  neces- 
sary, however,  to  stipulate  known  values  for  8w/8z  in  the  equation  of 
continuity  to  relax  these  conditions.  In  practice,  3w/3z  at  all  depths 
is  not  known,  but  models  with  "reasonable'1  values  are  considered  to 
test  the  restriction.  A less  drastic  assumption  is  neglecting  the  non- 
linear terms.  Checking  the  magnitude  of  these  te^ms  by  using  the  values 
of  velocity  that  resulted  from  the  computation  reveals  that  they  are  every- 
where less  by  order  of  magnitude'  than  the  terms  retained  in  the  equation 
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of  motion.  This,  in  part  is  due  to  the  low  resolution  resulting  from  the 
coarse  grid. 


The  conditions  of  no  velocities  at  the  lateral  boundaries  is  not 
correct  where  the  North  Atlantic  is  connected  with  other  bodies  of  water. 


The  results  at  these  places  should  be  considered  more  approximate  than 


those 

in  the  rest  of  the  ocean. 

In  spite  of  the 

obvious  crudity  of  the  model, 

it  is  felt  that  it 

can  serve  as  a 

useful  first  approximation  to  the  circulation. 

Table 

I:  Depth 
were 

s in  meters  of  levels  at  which 
computed . 

stream  functions  and  velocitie 

Level 

No. 

Depth 

Level  No. 

Depth 

1 

0 

17 

320 

2 

20 

18 

340 

3 

40 

19 

360 

4 

60 

20 

380 

5 

80 

21 

400 

6 

100 

22 

500 

7 

120 

23 

600 

8 

140 

24 

700 

9 

160 

25 

800 

10 

180 

26 

900 

11 

200 

27 

1000 

12 

220 

28 

1200 

13 

240 

29 

1400 

14 

260 

30 

1600 

15 

280 

31 

1800 

16 

300 

32 

2000 
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